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Ab  initio  quantum  chemical  simulations  were  coupled  with  variational  transition  state  theory  in  esti¬ 
mating  rate  constants  for  the  H  +  C3H3andH  +  C3H5  recombination  reactions.  The  energy  of  interaction 
between  the  H  atom  and  each  of  the  radicals  was  evaluated  at  the  CAS  +  1+2  level  for  the  range  of 
separations  and  relative  orientations  spanning  the  transition  state  region.  An  analytic  representation  of 
these  interaction  energies  was  then  implemented  in  variable  reaction  coordinate  transition  state  theory 
calculations  of  the  high  pressure  limit  recombination  rate  constant  for  temperatures  ranging  from  200  to 
2000  K.  For  the  propargyl  reaction,  the  overall  addition  rate  was  separated  into  contributions  correlating 
with  the  initial  formation  of  allene  and  propyne.  These  theoretical  results  were  compared  with  the  available 
experimental  data  as  well  as  with  corresponding  theoretical  estimates  for  the  H  +  C2H3  and  H  +  C2H5 
reactions.  The  H  +  propargyl  and  H  +  allyl  total  recombination  rates  were  remarkably  similar,  with  both 
being  greater  than  the  H  +  vinyl  and  H  +  ethyl  rates,  due  to  the  presence  of  twice  as  many  addition 
channels. 


Introduction 

It  has  been  widely  suggested  that  the  relatively 
high  concentrations  (in  the  combustion  environ¬ 
ment)  of  resonantly  stabilized  radicals,  such  as  the 
propargyl  radical,  are  the  result  of  their  compara¬ 
tively  low  rates  of  reaction  [1],  In  this  work  we  ex¬ 
amine  the  reaction  kinetics  for  the  recombination  of 
H  atoms  with  two  such  resonantly  stabilized  radicals; 
the  allyl  radical  (CH2CHCH2)  and  the  propargyl 
radical  (CHCCH2).  The  allyl  radical,  with  its  two 
degenerate  resonance  structures,  is  a  prototypical 
resonantly  stabilized  radical.  The  propargyl  radical 
also  has  two  dominant  resonance  structures,  but,  un¬ 
like  allyl,  the  two  structures  are  not  equivalent. 

There  is  considerable  recent  interest  in  the  rates 
of  propargyl  radical  reactions  due  to  the  postulated 
role  of  propargyl  radicals  in  the  soot  formation  pro¬ 
cess,  and  particularly  the  initial  formation  of  ben¬ 
zene  [1].  In  spite  of  this  interest,  the  available  ex¬ 
perimental  data  for  both  the  11+  propargyl  [2 — 4] 
and  the  II  +  allyl  [5]  reactions  are  quite  limited. 
The  present  coupling  of  high-level  quantum  chem¬ 
ical  simulations  with  sophisticated  transition  state 
theory  evaluations  should  provide  quantitatively 
meaningful  a  priori  estimates  for  the  kinetics  of 


these  two  reactions.  Comparison  with  related  theo¬ 
retical  results  for  the  H  +  vinyl  [6]  and  LI  +  ethyl 
[7]  reactions  allows  for  a  reasonably  direct  exami¬ 
nation  of  the  effects  of  resonance  stabilization  on  the 
kinetics. 

For  the  propargyl  radical  there  are  two  distinct 
addition  products,  allene  (CH2CCH2)  and  propyne 
(CH3CCH),  formed  by  attack  on  different  ends  of 
the  propargyl  radical.  For  allyl,  the  only  stable  ad¬ 
dition  product  is  propene  (CH3CHCH2),  which  can 
be  formed  by  attack  on  either  end.  Both  radicals  are 
planar,  and  attack  can  occur  from  either  above  or 
(equivalently)  below  the  plane  of  the  radical.  Thus, 
for  allyl  there  are  four  distinct  pathways,  all  of  which 
are  related  by  symmetry.  For  propargyl,  there  are 
also  four  distinct  pathways,  which  can  be  divided 
into  two  symmetry-related  pairs,  one  pair  leading  to 
allene  and  the  other  to  propyne. 

Here,  the  interaction  of  the  H  atom  with  each  of 
the  radicals  for  each  of  the  addition  paths  is  first 
characterized  with  high-level  ab  initio  quantum 
chemical  simulations.  The  resulting  analytical  poten¬ 
tials  are  then  implemented  in  variable  reaction  co¬ 
ordinate  transition  state  theory-based  estimates  for 
the  high-pressure  recombination  kinetics.  The  re¬ 
sults  of  these  calculations  are  directly  compared  with 
the  available  experimental  data. 
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ELEMENTARY  REACTION  KINETICS 
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Fig.  1.  Two-dimensional  cut  of  the  three-dimensional 
surface  for  the  H  +  C3H3  — »  C3H4  association  reaction. 
The  plotting  plane  is  perpendicular  to  the  plane  of  the 
C3H3  radical  and  bisects  the  HCH  angle.  Solid  contours 
are  positive,  dashed  contours  negative,  and  the  zero  energy 
contour  (defined  to  be  the  energy  of  the  H  +  C3H3)  is 
shown  with  a  heavy  solid  line.  The  contour  increment  is  2 
kcal/mol,  and  all  distances  are  shown  in  atomic  units  (one 
atomic  unit  =  0.52918  A). 

Potential  Surface  Calculations 

Methodology 

All  the  electronic  structure  results  reported  here 
are  from  multireference  configuration  interaction 
(MR-CI)  calculations  employing  orbitals  optimized 
using  the  complete  active  space,  self-consistent  field 
(CASSCF)  methodology.  In  these  calculations,  the 
CASSCF  reference  wave  function  consists  of  four 
active  orbitals  (the  hydrogen  Is  orbital  and  three  n 
orbitals  on  the  hydrocarbon  fragment)  and  four  ac¬ 
tive  electrons.  This  is  the  minimum  necessary  to  cor¬ 
rectly  describe  the  breaking  CH  bond  and  the  res¬ 
onance  of  the  hydrocarbon  radical.  The  effects  of 
higher-order  excitations  were  tested  using  a  multi¬ 
reference  Davidson  correction.  The  basis  set  used  in 
most  calculations  was  the  correlation-consistent,  po¬ 
larized  valence  double  zeta  (cc-pvdz)  basis  set  of 
Dunning  [8-10],  To  estimate  the  error  associated 
with  the  use  of  this  relatively  small  basis  set,  inter¬ 
nally  contracted,  MR-CI  calculations  were  carried 
out  along  the  association  paths  for  H  +  C3H3  using 
the  correlation-consistent,  polarized  valence  triple 
zeta  (cc-pvtz)  basis  set.  The  uncontracted  Cl  calcu¬ 
lations  were  carried  out  using  the  COLUMRUS 
package  of  codes  [11].  The  internally  contracted  Cl 
calculations  were  done  using  the  MOLPRO  package 


Results 

Characterization  of  the  potential  surface  for  the 
association  reactions  was  accomplished  by  mapping 
out  three-dimensional  grids,  in  polar  coordinates, 
corresponding  to  motion  of  the  hydrogen  atom  re¬ 
actant  around  the  rigid  hydrocarbon  radicals.  These 
surfaces  were  fit  using  techniques  described  previ¬ 
ously  in  calculations  on  II  +  C2H3  [6]  and  H  + 
C2H5  [7], 

A  two-dimensional  slice  of  the  three-dimensional 
H  +  C3H3  surface  is  shown  in  Fig.  1.  This  plot  il¬ 
lustrates  that  association  reactions  can  take  place 
through  four  barrierless  attractive  valleys:  two  sym¬ 
metry-related  paths  for  attack  on  the  CH  end  of 
C3H3,  and  two  for  attack  on  the  CH2  end.  Recall 
that  the  ground  state,  electronic  wave  function  of 
propargyl  is  an  unsymmetrical  resonance  between  a 
structure  having  two  carbon-carbon  double  bonds 
(with  the  radical  orbital  on  the  CH  end)  and  one 
having  one  carbon-carbon  single  bond  and  one  car¬ 
bon-carbon  triple  bond  (with  the  radical  on  the  CH2 
end).  The  C2v  equilibrium  structure  of  propargyl  im¬ 
plies  that  the  single  bond-triple  bond  structure  is 
the  dominant  resonance  structure.  Fig.  1  shows  that 
the  paths  for  attack  on  the  CH2  side  are  significantly 
more  attractive  and  wider  (in  this  in-plane  coordi¬ 
nate)  than  the  paths  for  attack  on  the  CH  side;  that 
is,  there  is,  not  surprisingly,  an  energetic  preference 
for  attack  of  the  dominant  resonance  structure.  In 
Fig.  2  we  show  two  perpendicular  slices  through  this 
potential  surface,  centered  on  each  of  the  terminal 
carbon  atoms.  From  this  plot  it  can  be  seen  that  the 
CH-side  paths  are  somewhat  wider  in  this  coordi¬ 
nate  than  the  CH2-side  paths.  This  is  caused  by  non- 
bonded  repulsions  between  the  incoming  H  atom 
and  the  two  CH  bonds  on  the  CH2  side. 

For  the  H  +  allyl  reaction,  the  calculations  show 
four  symmetry-related  pathways  for  barrierless  re¬ 
action  (the  hydrogen  can  attack  either  terminal  car¬ 
bon  from  either  above  or  below  the  plane  of  the 
radical).  The  calculations  also  predict  that  approach 
of  the  central  carbon  is  repulsive  (as  was  found  for 
H  +  C3H3). 

To  better  compare  the  relative  attractiveness  of 
the  association  paths,  we  optimized  minimum  en¬ 
ergy  paths  (MEP)  for  both  H  +  allyl  and  H  +  pro- 
pargyl.  In  these  calculations  the  C-H  distance  was 
kept  fixed  at  a  number  of  different  values,  and  the 
two  angular  coordinates  were  optimized  to  find  the 
orientation  with  the  lowest  energy  (again  the  ge¬ 
ometries  of  the  hydrocarbon  radical  fragments  were 
kept  fixed).  Plots  of  the  energy  along  these  MEPs 
are  shown  in  Fig.  3  along  with  comparable  calcula¬ 
tions  on  PI  +  C2H3  and  PI  +  C2H5.  From  this  plot 
it  can  be  seen  that  the  II  +  allyl  MEP  lies  between 
the  two  H  +  propargyl  MEPs.  All  of  the  H  +  res¬ 
onance  stabilized  radical  MEPs  were  predicted  to 
be  significantly  less  attractive  than  the  MEPs  for  the 
non-resonance  stabilized  radicals,  C2H3  and  C2H5. 


REACTION  OF  H  WITH  C3Hn  RADICALS 
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Fig.  2.  Two-dimensional  cuts  of  the  three-dimensional 
surface  for  the  H  +  C3H3  — »  C3H4  association  reaction. 
The  plotting  planes  are  perpendicular  to  the  CCC  axis.  The 
plotting  plane  for  the  upper  plot  contains  the  CH2  carbon, 
while  that  of  the  lower  plot  contains  the  CH  carbon.  The 
other  plotting  conventions  are  as  in  Fig.  1. 


The  MEPs  for  the  two  reactions,  H  +  C3H3  — » 
CH2CCH2  and  H  +  C3H3  — >  CH3CCH,  were  re¬ 
calculated  using  the  larger  cc-pvtz  basis  set.  The  cc- 
pvtz  MEPs  were  found  to  be  slightly  more  attractive 
than  the  cc-pvdz  MEPs.  The  cc-pvdz  potential  sur¬ 
face  was  corrected  for  this  deficiency  as  follows.  At 
each  point  the  distance  between  the  II  atom  and 
each  of  the  terminal  C  atoms  was  determined.  Two 
corrections  were  then  calculated,  one  using  the  dif¬ 
ference  between  the  cc-pvdz  and  cc-pvtz  MEPs  for 
addition  to  the  CH2  end  and  the  distance  between 
the  H  atom  and  the  CH2  carbon.  The  second  cor¬ 
rection  used  the  difference  between  the  two  MEPs 
for  addition  to  the  CH  end  and  the  distance  between 
the  H  atom  and  the  CH  carbon.  The  largest  of  these 


^C-H  W 


Fig.  3.  MEP  energies  for  the  reactions:  H  +  C3H5  — > 
C3H6  (solid),  H  +  C3H3  — >  CH2CCH2  (upper,  dashed),  H 
+  C3H3  — >  CH3CCH  (lower,  dashed),  H  +  C2H3  — >  C2H4 
(dash-dot),  and  H  +  C2H5  — >  C2H6  (dot). 

two  corrections  was  then  used  to  correct  the  cc-pvdz 
surface  at  that  point. 

Transition  State  Theory 

For  barrierless  reactions  such  as  the  radical-radi¬ 
cal  recombination  reactions  of  interest  here,  the  op¬ 
timal  form  for  the  transition  state  dividing  surface 
varies  considerably  with  energy  ( E )  and  with  total 
angular  momentum  (/).  At  low  E  and  low  J,  the  tran¬ 
sition  state  lies  at  large  interfragment  separations, 
and  the  optimal  reaction  coordinate  is  simply  the 
separation  between  the  centers  of  mass  of  the  two 
radicals.  At  higher  E  and/or  ] ,  the  transition  state 
moves  in  to  shorter  separations  where  the  optimal 
transition  state  dividing  surface  correlates,  at  least 
qualitatively,  with  a  fixed  separation  between  the 
centers  of  the  two  orbitals  involved  in  the  incipient 
bond. 

In  this  work,  both  the  definition  and  location  of 
the  transition  state  dividing  surface  were  variation- 
ally  optimized  within  the  set  of  surfaces  defined  by 
a  fixed  separation  between  the  H  atom  and  a  pivot 
point  connected  to  the  propargyl  or  allyl  radical. 
Choosing  the  pivot  point  as  the  center  of  mass  of  the 
radical  yielded  a  form  appropriate  for  low  E  and  J. 
Alternatively,  choosing  the  pivot  point  along  the  axis 
of  one  of  the  lone  pair  orbitals  of  the  radical  pro¬ 
vided  a  dividing  surface  that  was  appropriate  for  high 
E  and/ or  high  J.  Minimization  of  the  transition  state 
partition  function  with  respect  to  both  the  location 
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H  +  C3H3  -->  Products 


Fig.  4.  Plot  of  the  high-pressure  recombination  rate  con¬ 
stant  versus  temperature  for  H  +  C3H3,  including  the 
separate  contributions  corresponding  to  the  propyne  and 
allene  product  channels. 


of  the  radical  pivot  point  and  the  separation  between 
the  H  atom  and  the  pivot  point  then  provided  the 
present  variable  reaction  coordinate  transition  state 
theory  estimate  for  the  rate  constant. 

For  the  propargyl  radical,  there  are  four  distinct 
reaction  paths  for  the  addition,  corresponding  to  the 
addition  of  the  H  atom  from  either  above  or  below 
the  radical  plane  to  either  the  CIT  or  the  CH2  side 
of  the  radical.  The  presence  of  strongly  repulsive 
regions  separating  each  of  these  channels,  and  the 
extension  of  these  repulsive  regions  out  to  fairly 
large  separations,  suggests  that  it  is  appropriate  to 
optimize  separate  dividing  surfaces  for  each  of  these 
channels.  The  CH  (allene)  and  CH2  (propyne)  ad¬ 
dition  channels  were  distinguished  here  by  the  angle 
between  the  incoming  H  atom  and  the  CC  bond 
from  the  central  C  to  the  methylene  C.  Angles  of 
greater  than  90°  are  assumed  to  correlate  with  al¬ 
lene,  while  values  of  less  than  90°  are  assumed  to 
correlate  with  propyne. 

A  grid-based  procedure  was  employed  in  the  var¬ 
iational  minimization  at  the  E/J  resolved  level  for 
each  of  these  channels.  The  grid  of  radical  fixed- 
point  locations  was  defined  in  terms  of  the  vector 
from  the  fixed  point  to  the  C  atom  of  the  incipient 
bond.  This  vector  was  taken  to  lie  in  the  plane  per¬ 
pendicular  to  that  of  the  radical,  but  passing  through 
the  CCC  axis,  and  was  defined  in  terms  of  its  length 
d  and  angle  9  relative  to  the  CC  axis.  Sample  cal¬ 
culations  suggested  that  two  distinct  sets  of  ( d,6 ) 
should  be  considered.  The  first  set  was  centered 
about  (1.2,0)  and  effectively  yielded  the  center-of- 
mass  reaction  coordinate.  The  second  set  included 
d  values  ranging  from  0.0  to  1.25  and  angles  ranging 
from  about  70  to  110.  This  second  set  sampled  fixed 
points  in  the  neighborhood  of  the  radical  orbital. 
The  first  set  of  fixed-point  locations  was  coupled 


with  a  grid  of  H  atom  to  radical  fixed-point  separa¬ 
tions  R  ranging  from  about  5.0  to  about  3.25  A  with 
a  0.25  A  spacing.  The  second  set  of  fixed-point  lo¬ 
cations  was  coupled  with  a  grid  of  R  ranging  from 
about  3.0  to  1.75  A  again  with  a  0.25  A  spacing.  The 
VARIFLEX  software  package  [17]  was  used  in  eval¬ 
uating  the  number  of  available  states,  NEj  ( R,d,9 ) 
(i.e.,  the  E/J  resolved  transition  state  partition  func¬ 
tion)  for  each  point  in  this  grid,  and  in  performing 
the  Boltzmann  weighted  integral  over  E  and  J  of  the 
minimum  of  IVE j  on  the  grid  to  obtain  the  thermal 
rate  constant.  Interestingly,  the  optimized  dividing 
surfaces  were  again  found  to  correlate  closely  with 
the  contours  of  the  radical  orbitals  [6], 

F  or  the  allyl  radical,  there  were  again  four  distinct 
reaction  paths,  but  these  were  now  all  symmetrically 
related.  Each  of  these  channels  was  again  well  sepa¬ 
rated  by  repulsive  regions  of  the  potential;  thus,  op¬ 
timized  dividing  surfaces  were  considered  for  each 
channel  separately.  The  procedure  employed  in 
these  optimizations  is  closely  analogous  to  that  for 
the  propargyl  radical  and  is  not  reviewed  here. 

The  variable  reaction  coordinate  formalism  [18] 
employed  here  also  implements  an  assumed  sepa¬ 
ration  between  the  internal  vibrational  modes  of  the 
fragments,  labeled  the  conserved  modes,  and  the  re¬ 
maining  transitional  modes,  which  vary  in  character 
from  the  reactants  to  the  products.  Here,  the  prop¬ 
erties  of  the  conserved  modes  are  assumed  to  be 
invariant  in  the  process  from  reactants  to  the  tran¬ 
sition  state.  Then,  their  contribution  to  the  transition 
state  partition  function  cancels  the  corresponding 
contribution  to  the  reactant  partition  function.  In  re¬ 
ality,  both  the  vibrational  frequencies  and  the  opti¬ 
mal  geometry  will  vary  somewhat  as  the  fragments 
approach  one  another.  For  example,  the  H  atom(s) 
on  the  C  being  attacked  by  the  incoming  H  atom 
will  move  away  from  the  approaching  H  atom  by 
bending  out  of  the  plane  of  the  radical.  Concomitant 
with  this  change  in  bending  angle(s)  will  be  a  change 
in  its  vibrational  frequency  and  other  related  tor¬ 
sional  frequencies.  However,  such  variations  tend  to 
be  of  only  secondary  importance.  For  example,  for 
the  frontside  addition  channel  in  the  H  +  vinyl  re¬ 
action,  these  two  effects  yield  a  correction  of  only 
4%.  Thus  such  variations,  which  would  require  sub¬ 
stantial  additional  computational  effort,  were  ne¬ 
glected  here. 


Results  and  Discussion 

The  present  MR-CI-based  variable  reaction  co¬ 
ordinate  transition  state  theory  results  for  the  11  + 
propargyl  high-pressure  recombination  rate  constant 
are  plotted  in  Fig.  4.  The  ratio  of  the  propyne  to 
allene  formation  rate  constants  is  roughly  constant 
at  1.55  for  temperatures  up  to  1000  K  and  then  de¬ 
creases  gradually  with  increasing  temperature. 
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H  +  C  ,H.  ->  Products 


T  (K) 


Fig.  5.  Plot  of  the  high-pressure  recombination  rate  con¬ 
stant  versus  temperature  for  H  +  C3H5. 


H  +  Radical  — >  Products 


Fig.  6.  Plot  of  the  high-pressure  recombination  rate  con¬ 
stant  versus  temperature  for  H  +  C3H5  (solid),  H  +  C3H3 
(upper,  dashed),  H  +  C2H3  (dash-dot),  and  H  +  CoH5 
(dot).  The  two  lower  dashed  lines  are  for  the  propyne  and 
allene  (lowest)  channels  of  the  H  +  C3H3  recombination. 

reaching  a  value  of  1.45  at  2000  K.  This  larger  mag¬ 
nitude  for  the  propyne  formation  channel  is  readily 
explained  on  the  basis  of  the  greater  attractiveness 
along  the  MEP  (cf.  Fig.  3),  and  somewhat  larger 
range  of  angles  for  which  the  interaction  is  attractive 
for  this  channel  (cf.  Fig.  1). 

The  available  experimental  results,  which  are  in¬ 
directly  determined  from  kinetic  modeling  of  the  py¬ 
rolysis  of  toluene  [3],  and  the  oxidations  of  acetylene 
[2]  and  propargyl  [4],  also  are  plotted  in  Fig.  4.  The 
theoretical  predictions  are  in  quantitative  agreement 
with  the  most  recent  of  these  experimental  results, 
which  were  obtained  for  6.65  mbar  at  295  K.  Sample 
master  equation  simulations  of  the  pressure  depen¬ 
dence  suggest  that  the  discrepancy  of  a  factor  of  3 
to  10  with  the  other  experimental  results  may  be 
due,  at  least  in  part,  to  a  deviation  from  the  high- 
pressure  limit  in  these  experimental  studies. 


The  present  data  also  could  be  used  to  predict  CH 
dissociation  rates  in  allene  and  propyne,  for  which 
there  is  a  substantially  larger  body  of  experimental 
data  [19],  However,  as  has  been  well  described  by 
Kiefer  et  al.  [19],  the  proper  modeling  of  the  dis¬ 
sociation  process  requires  a  good  description  of  the 
effect  of  vibrational  anharmonicities  on  the  state 
density  of  the  complex,  which  is  beyond  the  scope 
of  the  present  work.  Our  own  sample  master  equa¬ 
tion  calculations  employing  the  present  VRC-TST 
model  for  the  transition  state  yielded  similar  results 
to  those  of  Ref.  [19],  Thus,  we  do  not  pursue  here 
a  comparison  with  dissociation  rates. 

Related  calculations  of  the  H  +  C3H3  recombi¬ 
nation  rate,  employing  a  potential  which  does  not 
include  the  cc-pvtz  basis  set  corrections,  yielded  re¬ 
sults  that  were  only  9%  to  19%  lower.  Similarly,  the 
neglect  of  the  Davidson  correction  yielded  results 
that  were  only  8%  to  20%  lower.  The  smallness  of 
these  two  corrections  suggests  that  any  inaccuracies 
in  the  underlying  potential  are  likely  to  be  of  only 
minor  significance. 

The  present  results  for  the  H  +  allyl  recombina¬ 
tion  rate  are  plotted  in  Fig.  5.  The  theoretical  pre¬ 
dictions  are  in  quantitative  agreement  with  the  fairly 
direct  experimental  measurement  of  Hanning-Lee 
and  Pilling  [5],  The  predicted  rise  by  a  factor  of  1.54 
from  200  to  2000  K  is  identical  to  that  predicted  for 
the  H  +  propargyl  reaction.  A  cc-pvtz  calculation  of 
the  MEP  was  not  performed  for  this  reaction.  How¬ 
ever,  the  near  constancy  of  the  effect  of  this  correc¬ 
tion  on  the  rates  for  a  variety  of  reactions  led  us  to 
simply  incorporate  a  15%  increase  to  the  rates.  The 
effect  of  the  Davidson  correction  is  also  expected  to 
be  similar  to  that  for  the  H  +  propargyl  reaction. 

A  plot  of  the  VRC-TST  predictions  for  the  high- 
pressure  recombination  rate  constants  for  the  H  + 
propargyl,  allyl,  vinyl  and  ethyl  reactions  is  provided 
in  Fig.  6.  The  rate  constant  for  the  H  +  C3H5  re¬ 
combination  is  remarkably  similar  to  the  total  re¬ 
combination  rate  constant  for  H  +  C3H3.  Interest¬ 
ingly,  the  H  +  C3HS  MEP  lies  roughly  midway 
between  the  MEPs  for  the  two  II  +  C3H3  channels. 
As  a  result,  while  the  addition  to  the  CH2  side  of 
propargyl  occurs  at  a  greater  rate  than  does  the  ad¬ 
dition  to  one  side  of  allyl,  the  addition  to  the  CH 
side  of  propargyl  occurs  at  a  lower  rate.  The  sum  of 
the  rates  for  the  two  propargyl  channels  then  adds 
to  be  essentially  identical  to  the  total  allyl  rate. 

The  H-atom  recombination  rates  for  the  two  res¬ 
onantly  stabilized  radicals  were  actually  predicted  to 
be  greater  than  those  for  the  vinyl  and  ethyl  radicals. 
At  first  glance  this  seems  to  contradict  the  increased 
attractiveness  for  the  MEPs  in  the  vinyl  and  ethyl 
radicals  (cf.  Fig.  3).  However,  this  result  is  simply  an 
indication  that  there  are  essentially  four  separate  ad¬ 
dition  channels  for  each  of  the  resonantly  stabilized 
channels,  while  there  are  only  two  (frontside  and 
backside  additions)  for  the  vinyl  and  ethyl  radicals. 
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TABLE  1 

Modified  Arrhenius  parameters  for  H  +  radical  high 
pressure  recombination  rates 


Reaction 

A  (10  11 

cm3  s_1) 

n 

Ea  (K) 

H  +  C3H5  — >  products 

5.804 

0.2363 

-25.90 

H  +  C3H3  — >  products 

4.520 

0.2711 

-42.24 

H  +  C3H3  — >  propyne 

3.710 

0.2298 

-24.26 

H  +  C3H3  — >  allene 

1.1.39 

0.3320 

-  69.06 

H  +  C2H3  — >  products 

4.279 

0.2724 

-37.56 

H  +  C2H5  — >  products 

10.38 

0.1017 

15.84 

After  division  by  two,  the  rates  for  the  allyl  and  pro- 
pargyl  radicals  were,  as  expected,  significantly 
smaller  than  those  for  the  vinyl  and  ethyl  radicals. 

Summary 

The  H  +  allyl  and  II  +  propargyl  recombination 
reactions  were  studied  with  variable  reaction  coor¬ 
dinate  transition  state  theory  employing  interaction 
potentials  obtained  from  multireference  configura¬ 
tion  interaction  calculations.  The  predictions  for  the 
high  pressure  limit  recombination  rates  were  found 
to  be  essentially  identical  for  these  two  reactions. 
The  results  for  the  H  +  allyl  reaction  were  in  good 
agreement  with  the  only  experimental  study,  while 
those  for  the  H  +  propargyl  reaction  agreed  only 
with  those  from  the  most  recent  experimental  study, 
being  considerably  greater  than  other  earlier  studies. 
The  addition  to  the  CH2  side  of  the  propargyl  radical 
was  predicted  to  occur  with  a  rate  that  is  1.5  times 
that  for  the  CH  side. 

Table  1  provides  the  parameters  A,  n,  and  Ea  from 
modified  Arrhenius  fits  [k  =  A  T"  exp (-EJT)]  to 
the  high-pressure  recombination  rates  for  H  +  allyl, 
H  +  propargyl,  H  +  vinyl,  and  H  +  ethyl.  Contrary 
to  initial  expectations,  the  high-pressure  recombi¬ 
nation  rates  for  the  two  resonantly  stabilized  radicals 
studied  here  were  actually  greater  than  those  for  the 
H  +  vinyl  and  H  +  ethyl  reactions,  due  to  the  pres¬ 
ence  of  twice  as  many  addition  channels  for  the  res¬ 
onantly  stabilized  radicals.  Similar  findings  might  be 
expected  for  the  relative  recombination  rates  with 
other  addition  partners. 

The  observation  of  elevated  concentrations  of 
these  radicals  thus  appears  to  be  simply  the  result  of 
greater  rates  of  formation;  that  is,  their  stabilization 
energy  should  correlate  with  decreased  production 
thresholds  relative  to  other  related  channels.  At 
lower  pressures  these  decreased  production  thresh¬ 
olds  might  also  affect  the  ordering  of  various  effec¬ 
tive  recombination  rates.  In  particular,  the  redisso¬ 
ciation  of  the  addition  products  back  to  reactants 
would  be  expected  to  be  greater  for  the  resonantly 


stabilized  radicals  due  to  their  lower  binding  energy, 
and  thus  the  effective  loss  rates  would  be  smaller. 
Of  course,  such  a  pressure-dependent  effect  would 
also  depend  quite  strongly  on  the  lifetime  (and,  cor¬ 
respondingly,  the  size)  of  the  complex. 
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COMMENTS 


Alexander  Fridman,  University  of  Illinois  at  Chicago, 
USA.  The  theoretical  modeling  presented  shows  an  essen¬ 
tial  positive  activation  energy  of  the  hydrogen-propargyl  H 
+  C3H3  recombination.  Recent  experiments  give  no  acti¬ 
vation  energy,  actually  no  temperature  dependence  at  all 
for  the  reaction.  Generally  similar  recombination  processes 
have  no  essential  positive  activation  energies.  Probably, 
such  a  discrepancy  is  due  to  approximate  nature  of  the 
approach. 

Author’s  Reply.  Current  experiments  are  probably  not 
sufficiently  accurate  to  resolve  the  question  as  to  whether 
or  not  there  is  a  small  positive  activation  energy  for  this 
reaction  in  the  high  pressure  limit  as  predicted  by  our  cal¬ 
culations.  Our  approach  does  give  a  negative  temperature 
dependence  for  other  reactions  where  the  experimental 
evidence  for  a  negative  temperature  is  stronger.  See,  for 
example,  CH3  +  CH3  [1]. 
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1.  Klippenstein,  S.  J.,  and  Harding,  L.  B .,/.  Phys.  Chem. 
103:9388-9398  (1999). 

• 

Jurgen  Troe,  University  of  Gottingen,  Germany.  The 
transition  state  theory,  which  you  use,  is  an  adiabatic  the¬ 
ory,  that  is,  it  corresponds  to  the  limiting  case  of  adiabatic- 
capture  dynamics.  We  have  shown,  however,  that  high 
pressure  combination  reaction  of  hydrogen  atoms  with  rad¬ 
icals  (or  other  open  shell  species)  are  showing  partly  non- 
adiabatic  dynamics  [1,2],  I,  therefore,  estimate  that  your 
results  for  200  K  have  to  be  reduced  by  about  20%  and  for 
300  K  by  about  10%  due  to  these  rovibrationally  non-adi- 
abatic  dynamical  effects. 
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2.  Maergoiz,  A.  I.,  Nikitin,  E.  E.,  Troe,  J.,  and  Ushakov, 
V.  G .,/.  Chem.  Phys.  108:5265,  9987  (1998). 

Author’s  Reply.  A  statistical,  rather  than  adiabatic,  theory 
was  used  in  this  work.  For  a  reaction  coordinate  defined 
by  the  separation  of  the  centers-of-mass  of  the  two  reacting 
fragments,  variational  statistical  theory  can  effectively  be 
reduced  to  an  adiabatic  channel  theory.  However,  for  gen¬ 
eralized  reaction  coordinates,  such  as  those  used  here, 
there  are  kinetic  couplings  between  the  reaction  coordinate 


and  the  remaining  coordinates.  These  couplings  introduce 
significant  conceptual  differences  between  the  statistical 
and  adiabatic  channel  frameworks.  The  correct  treatment 
of  these  couplings  is  an  important  component  of  the  cur¬ 
rent  statistical  theory.  Importantly,  it  should  be  possible  to 
recover  the  effect  of  the  non-adiabaticities  in  the  center- 
of-mass  frame  via  the  appropriate  choice  of  generalized 
reaction  coordinate  within  a  variational  statistical  theory 
framework.  Thus,  one  cannot  use  the  comparison  of  tra¬ 
jectory  simulations  with  adiabatic  capture  theory  to  suggest 
that  the  present  results  should  be  reduced  by  10%  to  20%. 


P.  R.  Westmoreland,  University  of  Massachusetts,  Am¬ 
herst,  USA.  We  have  proposed  that  propargyl  +  propargyl 
combination  should  favor  combination  at  the  CH2  ends,  as 
you  find  for  the  site  of  H  +  propargyl  |  lj.  Our  simplistic 
model  was  that  if  90%  of  the  electron  density  was  on  the 
CH2  carbon  and  10%  on  the  CH  carbon  (calculated),  the 
probabilities  of  l,5-hexadiyne:hexa-4,5-dienyne:l,2,4,5- 
hexatetraene  would  be  0.81:0.18:0.01.  [Note  that  the  hex- 
adiyne  and  hexatetraene  can  interconvert  by  a  Cope  rear¬ 
rangement.]  Wouldn't  your  results  support  favoring  the 
hexadiyne? 
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Author’s  Reply.  Our  calculations  predict  the  CH2  end  of 
propargyl  is  more  reactive  (toward  hydrogen  atoms)  than 
the  CH  end  by  a  ratio  of  0. 6:0.4.  If  we  assume  that  this 
difference  in  the  reactivity  of  the  two  ends  of  propargyl 
carries  over  directly  from  H  +  C3H3  to  C3H3  +  C3H3 
then  we  would  predict  the  three  products  you  mention,  1,5 
hexadiyne,  hexa-4,5-dienyne  and  1,2, 4, 5  hexatetraene,  to 
be  formed  in  the  ratios  0.36:0.48:0.16  respectively.  How¬ 
ever,  it  is  generally  found  that  less  reactive  species  will  react 
more  selectively  than  more  reactive  species.  Since  propar¬ 
gyl  is  certainly  less  reactive  than  a  hydrogen  atom,  we 
would  expect  the  above  prediction  of  36%  1,5  hexadiyne 
to  be  a  lower  bound.  In  fact,  recent  measurements  con¬ 
cluded  that  the  yield  of  1,5  hexadiyne  from  this  reaction  is 
60%  [1], 
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